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Introduction. A special case of synthesizing the same 
electromechanical control system by the Pontryagin maximum 
principle and by the synergetic synthesis method is considered. 
The task was to solve the synthesis problem of the time 
optimal electromechanical position control system; herewith 
the travel resistance modulus linearly depended on the output 
coordinate of the system. This approach to the selection of the 
synthesis problem was because the synthesis of time optimal 
systems is one of the most widespread problems, and it is 
solved by increasing the efficiency of the existing control 
systems. 

Materials and Methods. Synthesis of the time optimal linear 
control system based on the maximum principle is a widely 
accepted problem in the modem control theory. However, the 
procedure of synergistic synthesis does not have such 
formalization. This being the case, the paper suggests an 
approach that brings together these two methods, which, in our 
opinion, will increase the efficiency of the synergistic 
synthesis method through adding some features of the 
synthesis methodology for optimal systems. 

Research Results. The paper formulates two key concepts. The 
first one is as follows: the application of the maximum 
principle for an object of the DC motor class when 
synthesizing the positioning algorithm under the conditions of 
linear loading functionally dependent on the engine rotation 
angle allows the time optimal system to be optimized. The 
second concept states that synthesis of a control system based 
on the synergistic approach enables to obtain a system close to 
optimal (quasioptimal), but after modifying the synergetic 
synthesis method itself. A hypothesis is formulated on the 
connection between the introduced 


possible (when 


implementing the procedure of state space extension in the 
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Beedenue. Paccmotpeu yacTHbIi ciryyalt CHHTe3a ODHOM UH Tor 
%Ke 9JICKTPOMexXaHH4ecKOH CHCTeMbI ylpaBJIeHHA MeTOJOM 
MakcuMyMa IlouTpsrHHa HM MeTOJOM cCHHepreTH4ecKoro 
cuute3a. B xauyecTBe 3aqauH Obila ompeyemena 3ayqaya 
CHHTe3a ONTHMaJIbHOL 10 OxicTpoelicTBHIO 
3JICEKTPOMeXaHH4eCKOH CHCTeMbI MO3HUMOHHpOBaHHA, pH 
39TOM MOMCHT COMPOTHBJICHHA JBYKCHHIO JIMHeMHO 3aBHCell 
OT BbIXOAHOM KOOPAHHATEI CHCTeMBI. DTOT MOAXOA K BbIOOpy 
3aqauH cHHTe3a Obi oOycNOBIeH TeM, 4TO CHHTEe3 
ONTHMAJIbHBIX T10 ObICTpOelCTBHIO CHCTeM ABJIAeTCA OMHOM 
43 CaMBIX WIMpOKO paciipocrpaHeHHbIx 3ajau, KOTOpad 
pelaeTca pH MOBbILIeHHH 9Pd@eKTHBHOCTH AelcTByIOWHXx 
CHCTeM ynpaBleHHa. 
Mamepuanst u  memodol. CuatTe3 OnNTHMaJIbHOM m0 
ObIcTpoyelcTBHIO IMHeHHOM CHCTeMbI ypaBeHHA Ha OCHOBe 
IIpHHHlla MakKCHMyMa — UIMpoKo pacipoctpaHeHHad 3ayaya 
B COBpeMeHHOM TeopHH ynpaBsieHua. OHako lpoleyypa 
CHHepreTH4eckoro CHHTe3a Tako (opMamM3aluH He MMeeT. 
Ucxoya “3 STOrO, B CTaTbe NpeWoxKeH NOAXoOL, cOmMKaIouIMi 
9TH JBa MeTOa, KOTOPbIM, 10 MHCHHIO aBTOPOB, NO3BOJIMT 
TIOBLICHTh 230p(PeKTHBHOCTb MeTOJa cHHepreTH4ecKoro 
CHHTe3a, 00aBMB B HerO HeKOTOpble ocoOeHHOCTH 
MeCTOJOJIO“MM CHHTe3a ONTHMAJIbHBIX CHCTeM. 

Pe3yiemamei uccredoeanun. B padote copMyMpoBaHbl Ba 
OCHOBHBIX Hay4HbIX TWoOKeHHA. Ileppoe - tpumMeHeHve 
TIpHHYWia MaKCHMyMa JVId OObeKTa kKylacca JBHTaTelib 
MOCTOAHHOTO TOKa mpu CHHTe3e 


TIOSHUMOHHpPOBaHHA = B YCJIOBHAX JIMHeHHOH 


amIropHt™ma 
Harpy3KH, 
(YHKUMOHAIbHO 3aBMcAlel OT yrla WoBopota ABHraTess, 
TIO3BOJIT€T ONTHMH3HPOBaTb CHCTeMy 0 ObICTpoyelicTBHIO. 
Bropoe - CHHTe3 CHCTeMBI ylIpaBeHHt Ha OCHOBe 
CHHepreTH4ecKkoro MOAXOa WO3BONAeT MOWYYHTh CHCTeMy, 
OnM3KYIO K OMNTHMaIbHOM (KBa3HONTHMAaJIbHYIO), HO yrke 
ocle MOMpuKaWHH caMoro MeTOa cHHepreTHuecKoro 
cuHTe3a. CcopmysMpoBaHa THoTe3a O BO3MO%KHOM CBA3H 
Mey peamm3alHH  mpolexyppr 
pacuimpenusa COCTOAHHA = B  MeTOze 


BBOJMMBIMH,  TIpvi 
IIpoctpaHcTBa 
CHHepreTw4ecKOrO CHHTe3a, MOCTOAHHBIMH BpeMeHH Cc 
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synergetic synthesis method) time constants with the optimal onpeyeaeMbIM B  MeTO€ MaKCHMyMa _ OIITHMaJIbHbIM 
switching time of control defined in the maximum method. BPCMCHEM TEP CKINOYCHHA YIPAaBICHH A. 

Discussion and Conclusions. The synthesis through the Coc prcetie i 3aicHo tenis: BrioureH cuales yupaEneHia 
maximum control technique and the ADAR method is MeTOJOM MaKCHMyMa 4H MeTOqOM AnasINTH4ecKOrO 
performed. In virtue of the comparison of efficiency of these Kouctpyupopanua ArperuposaHHbix Perynatopos AKAP. Ha 
methods, a hypothesis is put forward on the possible OCHOBaHHM cpaBHeHua 9:)eKTHBHOCTH IIpPHMeHeHHA MeTOOB 


compatibility of the studied methods. BbIIBHraeTcA THMOTe3a O BO3MO2%XKHOM COBMeCTHMOCTH 


MCCIICHYECMBIX MCTOJIOB. 


Keywords: maximum principle, optimal control, operation KuroyveBbie ClOBa: TpHHIMI MaKCHMyMa, OITHMaJibHoe 
speed, control algorithm, synergetic synthesis, ADAR method. ympaBienve, ObICTpoyelicTBHe, aJIOpHTM ylipaBeHus, 
ciHepretTnyeckuii cuHtTe3, MeTtoa AKAP. 
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Introduction. In the 20th century, the engineering requirements, in particular, in space engineering, put 
forward a range of problems for which a new theory was developed — the theory of optimal control [1]. One of the main 
techniques for the synthesis of optimal control systems is the maximum principle developed by Soviet mathematician L. 
S. Pontryagin and his disciples in the fifties-sixties of the 20th century [2]. The application of this principle is based on 
the formalization of the synthesis problem with the transition to the form of the Mayer problem and the subsequent 
solution to systems of linear or, in some special cases, nonlinear differential equations [1]. It is worth noting that the 
need to solve systems of differential equations, and in the nonlinear case this is not always possible, is in many ways a 
limitation of both the maximum principle itself and the whole concept of synthesis of optimal control systems. 

In recent decades, a new approach to the synthesis of control systems has been widely adopted. It relies on the 
synergetic concept of the analysis and synthesis of systems. The technique used, the author of which is A. A. 
Kolesnikov, is called the method of analytical design of aggregated regulators (ADAR) [4-5]. The proposed approach is 
based on the concept of synthesis of nonlinear feedbacks. They provide the asymptotic stability of the control system 
with respect to the required motion of the attractor in the state space [4]. This method differs from the methods of 
synthesis of optimal control systems in the absence of both the optimization criterion of the control system and the 
statement of the synthesized control optimality. From the point of view of the implementation of the synthesis 
procedure, the ADAR method has an undeniable advantage over the synthesis methods of optimal systems, which is 
expressed in the absence of restrictions on the nonlinearity of the system of differential equations [4-5]. 

From a practical standpoint, synthesis of the process or object management should be able to answer the 
question if there is another control that has the property of superiority over all others. Thus, the modern mass production 
constantly requires efficiency improvement to ensure market competition. One of the most popular ways to increase this 
efficiency is to optimize management processes. This approach is applicable both to the systems of automated assembly 
of equipment [6—7] and to the metal-cutting systems, in particular, to drilling control systems in metal-cutting machines 
[8-15]. The idea of combining these approaches to synthesis of the control systems has, from the authors’ the point of 
view, an undoubted practical effect. In science terms, it is important to combine the ADAR method advantages, which 
are expressed in the possibility of considering synthesis of the complex nonlinear process dynamics in mechanical 
engineering [16-19], with a neat and definite formalization of the synthesis problem formulation and assessment of its 
achievement in the maximum principle [1—2] . 


1. Synthesis of basic mathematical model and formulation of research problem 

In the modern economy, the direction of the time optimal system synthesis, which enables to obtain a 
significant increase in the global system efficiency [7-8], has become a frequent practice. Thus, under solving the 
problems of automating the assembly processes of various equipment, the task of attaching different types of parts to 
each other [6—7] often arises. Here, the economic efficiency of the entire production process depends on the speed of 
this operation. The same situation is observed in metal cutting systems on metal-cutting machines [9-19], in which the 
faster the machining process is, the lower the costs of the entire production process. Based on the reasoning, under the 
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assumption that it is necessary to synthesize a time optimal control system, which, considering the similarity of the 
processes of automated fastening of parts and drilling deep holes, can be illustrated by the following diagram (Fig. 1). 


C 
DC motor 


Fig. 1. Illustration diagram 


Fig. 1 shows a DC motor performing either an equipment fastening operation or a feed in a drilling unit. In the 
first case, C is the combined load-deflection characteristic of the material that is compressed under twisting. In the 
second case, C characterizes a linear increase in the feed resistance under the accumulation of chips in the chip removal 
channels of the drill. 

The electric drive, which ensures the conversion of electrical energy to mechanical energy of motion of the 
system actuators, is the basis for the support subsystem of the control system in both cases considered. Let us assume 
that in the present case, such a transformation is carried out by a DC motor with collector control, which is described 
through the following system of equations [20]: 


where U is voltage applied to the engine manifold; i is current consumed by the motor; R, L are parameters of the 
electrical part of the engine; J is parameter characterizing the inertial properties of the engine rotor, the reduced 
inertial moment of all rotating masses; w is engine rotor speed; Mis external moment of resistance; Cm, C. are 
mechanical and electric engine constants. In this case, the moment is presented as a function of the angle of rotation of 
the engine rotor, that is M, = Ca . With this in mind (1), we rewrite as: 


Uetese we 
dt 


pea ace ¢ (2) 
dt 

da 

—=0 

dt 


We express the current value in the second equation, apply the obtained value in the first equation, and get: 
1 LI d’o RJ dw LC RC 
U O= st + o+ OL 
c. Cc. Cn dt c. Cy dt C. Cn C. Cn a (3 ) 
da 
“=o 
dt 
We solve the first equation with respect to the highest derivative and rewrite it with the second one, and then 


we receive: 

















da 
ds 
do Rdo «4c C.. RC... ¢,.. (4) 
= ( +—)o a+—“U 
dt Ldt LJ J LJ LJ 
d : R R 
Let us puta =x,,@= 4 =x,; as constants, we introduce — = Oy (—2 c+ = = days oe a,,,—>=b. 
dt L LI J LI LJ 


Then the system (4) takes the form: 


Vestnik of Don State Technical University. 2018. Vol. 18, no. 4, pp. 438-448. ISSN 1992-5980 eISSN 1992-6006 
Becmuuk Jloncxozo zocyOapcméennozo mexnuyeckozo yHueepcumema. 2018. T. 18, Ne 4. C. 438-448. ISSN 1992-5980 eISSN 1992-6006 








ca 
i 
Xy 
— =x ‘ 5 
“i 3 (5) 
IX. 
a = A,X, — Ay X, — AX, +bU 
or in matrix-vector form: 
x= Ax+Bu, (6) 
0 1 0 0 Xj 0 
where A=| 0 0 1 |, B=|0}, x=| x, |,u=] 0 
743, —43, 33 b x; U 


We take the DC motor whose parameters provide the following constant values: a3;=4.65, a3.=4.6, a33=2, as a 
drive ensuring the system motion. 
0 1 0 
A matrix, with this in mind, will take the form: A=| 0 0 1 |. From knowing A matrix, we find 
4,65 -4,6 -2 


eigenvalues of A matrix: 2, = —1.2666 + 0.0000i, 1, = —0.3667 + 1.88061, 13 =—0.3667 — 1.88061. 

As is clear from the obtained eigenvalues of the matrix, the control system is asymptotically stable according 
to Lyapunov [21]. Moreover, we can argue about the oscillatory nature of the processes proceeding in the system, since 
the eigenvalues contain not only negative real parts, but nonzero imaginary parts. 

2. Synthesis of control by Pontryagin’s maximum principle 
First, we formulate the problem of the optimal control synthesis in the following form: 

- using the maximum principle for the object described by the System (5), determine the optimal equation 
algorithm that ensures the transfer of the object from the initial state x,(0)=50, x,(0)=0, x;(0)=0 to the final state 
x)(T)=0, x2(T)=0, x3(T)=0 for T minimum time. Herewith, it is necessary to determine transition count, switching 
torque, and to construct curves of u(t) control and x/(t), x2(t), x3(t) coordinates. 








lulSU max Yestriction is imposed on the control action. The system parameters are as follows: 
a,, =1,a,, =12,a,, = 20,b=1, Umax=440 V. 


Solution: We write the Hamiltonian: 





H= WSs: 7) 


where y, and f; are coordinates of wy and f vectors. In addition to the system (5), we construct a system of equations 


for y, auxiliary variables, where i = 1,2 , using the following relation: 


dy, Of, 
= ot fai: 
dt DY: ox ‘ 


j=l i 





Or open : 


dw, of of, a; 
dt E We Be ae 




















WM eh 1 + of, 2 + Is 3 (8) 
dt OX, Ox, Ox, 
dy, oh a, Hs 
+ + 
dt E ge a 
Considering (5), the equation system (8) takes the form: 
oe =4,65¥, 
dt 

TO =, 44,69, (9) 

a =, 42, 


dt 
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d¥, 








=4,65, 
+ 
— =-Y, +46, (9) 
— =-Y,+2¥, 
t 
Eigenvalue matrix for (9) case: 
1.2666 0 0 
a =) © 0.3667 +1.8806i 0 
0 0 0.3667 + 1.8806: 
The solution to the system (9) for the diagonalized case of AY? matrix takes the form: 
WP. = Ce! 
janice | 
Y, = ae : (10) 


_ (0.3667-1.8806i)¢ 
Y, =Cye 


Considering V¥ matrix, the solution in the initial basis will be: 
W, =-0.9474C,e' + (0.5812 — 0.33471) C, e077 18" + (0.5812 + 0.33471), e077 0" 
oe = -] : ] 893C,e'7°" + 0.6883, 0300771 880601 = 0.68e3C ee 
W, =—0.2581C,e'7" + (0.1812 + 0.20861)C, 60777! 8 +. (0.1812 — 0.20861) C, 6077 1 80" 
where C,, C2, C3 are integration constants. 
The general expression describing the Hamilton function: 
H=¥f+¥, 4+ (11) 
In the expression (11), an important — from the point of view of the method of synthesis — role is played by the 
member, which includes the control: 
H’ =U =(-C,0.2581e'*™ + C, (0.1812 + 0.20867) + C, (0.1812 — 0.20861) e077 "bu 
In order for H Hamiltonian calculated by the formula (11) to take the maximum positive value, H’ term must 
be always positive and the greatest. For this, the optimal control algorithm should be u(t)=OU nq,, where: 
o = sign(—0.2581C,e'7 + (0.1812 + 0.20861) e077!" + (0.1812 — 0.20867) E077 80) 


u(t) optimal control is a piecewise constant function, taking +U,,,, values, and it has no more than two 





intervals of constancy, since the nonlinear function 

—0.2581C,e'7°" + (0.1812 + 0.20867) e071" +. (0.1812 — 0.20867) e007 189" 
changes the sign no more than once. In this case, the possible sign change occurs from plus to minus, that is, to fulfill 
the maximum principle, it is required to first apply U= +U,,,x, to the engine, and then U= -U,,,,. Let us verify these 
arguments through constructing ‘¥, obtained functional dependence in the Matlab package (Fig. 2). 








Fig. 2. Graph of ‘V, variance 
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As can be seen from Fig. 2, ¥, does change the sign from (+) to (—) only once. With this in mind, we define an 


optimal equation algorithm that ensures the transfer of an object from x,(0)=0, x2(0)=0, x3(0)=0 initial state to 
x,(T)=50, x2(T)=0, x3(T)=0 final state for 7 minimum time. Here we note that under modern conditions, there is no 
need to obtain an analytical solution to the original system of equations. Using the available capabilities of modern 
software packages like Matlab, we straightforward and clearly can obtain a numerical solution to the case under 
consideration. 

The numerical simulation results for the calculated time optimal control algorithm are presented in Fig. 3. The 
model parameters are selected in such a way that it fulfills the boundary condition required by the right-hand border. 


Xi, rad 





0 05 1 15 tsec 








0 05 1 1.5 
b) 


Fig. 3. Graphs of output coordinate (a) and control switching 
at ¢, =0.5(b) 

Fig. 3 shows that the control switching time from (+U max) to (-Umax) is selected as t, = 0,5 seconds, and the 
total control time is T=1.3 seconds. In this case, the control system comes to the required output level, that is, the 
right-hand boundary of the boundary conditions is reached. 

3. Synthesis of control system through ADAR method 
To synthesize a control system using the ADAR method, it is possible not to make a transition to the abstract 


case of the state space (see (5)), but it is easy to use the original system (2). However, it needs to be translated to the 
following form: 








dx, _ 
Lo 
X, 
2 
hi = Ay) X, + Ay3X3 2 (12) 
dx. 
ae = Gy Xy — Ay, X, +bU 
t 
. . Ch, C C, R 1 : 
where Q@=xX,,@=X,,i=x,are variables; = 43, a>), a 33, b are introduced as constants. 
i} 2 3 J 23, J 21 L 32, L 33 L 


Considering the previously defined values, we receive: a,, =1.515,a,, =2.325,a,. =1.5,a,,=2,b5=1 and the system 
(12) takes the form: 


a 

fe 

= 72.325%, +1.515x, , ae 
dx, _ —1.5x, -2x,+U 

dt 


With this in mind, A matrix takes the form: 
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0 1 0 
A=|-2.325 0 1.515}. 
0 -15 -2 
From knowing A matrix, we find A eigenvalues of A matrix: A1 = —1.2674 + 0.00001; A2 = —0.3663 + 1.88011; 
13 = —0.3663 - 1.8801i. As is clear from these results, it is referred to the control system case considered in the 


synthesis through the maximum principle. 

The coordinate characterizing the angle of rotor rotation (x ;) is the output coordinate of the system. Therefore, 
in order to form requirements for the desired system behavior in the state space, we introduce a macrovariable of 
Y =x,—X), = Vorder, where x,, is the specified target value of x, coordinate. The velocity value of the angle change 


must ensure that the angle of rotor rotation tends to the value specified above. To this end, we introduce another 
macrovariable — ‘Y= x,—6,(x,)=>0, where,(x,) is a certain function describing x, => ,(x,) tendency in the 
stationary state. Further, we introduce a macrovariable of order in Y =x, -6,(%,,x,)=>0 coordinate, where 
;(x,,xX,) is a certain function that describes the relationship between the coordinates in the stationary state of the 
system and, therefore x, > ,(%,,x,). 

For newly introduced ¥,, ‘¥,, ‘¥, macrovariables, we require the asymptotically stable law of change, that is, 


execution of the following system of differential equations: 











7 +¥,=0 

t 
d¥, 

T, , +¥,=0, (14) 
Y 

74 pete! 


where 7,, T,, T, are time constants that determine the vanishing rate of the introduced macrovariables, or, in other 
words, the shrinking rate of an arbitrary path of the system in the state space to the desired attractor. 

The system of equations introduced in the expression (14) has expanded the system state space from the 3rd to the 
6th, which is one of the main provisions of the ADAR method. Further synthesis of the control system is reduced to a 
stage-by-stage process of decomposition (compression) of this space to the initial level [1-2]. We start the 
decomposition process with the substitution of the system (14) to the first equation instead of ‘¥; macrovariable of its 
value in the coordinates of the controlled process. As a result, we obtain: 


dx, 
A to) = 9, (15) 
Considering (13), the equation (15) is written as: 
Tx, +(x, —Xo,) = 90. (16) 


Considering the next step of the synthesis algorithm, because of which it was assumed that x > o,(x,), we find 


, (x, ) value as: 


(% =o) 
$,(x,) =-— (17) 
T, 
With this context, ‘¥, macrovariable is as follows: 
W, =x, -0,(¢) =, + a (18) 


1 
Validity of the expression (17), in terms of the control objectives, is confirmed by the fact that the value in the steady 


state is as follows: x, > 6,(x,) > 0. Then, considering (18), the second equation of the system (14) will take the 


following form: 


dx do,(x,) 
Te a) =O (19) 
Or, considering ,(x,): 
: dx, re dx, gage pd _ 
dt T dt ~ T, 





0. (20) 
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At the next step of the synthesis algorithm, using the same reasoning as before, and the fact that we 


determine x, > ,(x,,x,) from (20), considering (13), ,(x,,x,) value in coordinates of the controlled process: 

T,+T, — 

( 2+ We 0,66 Xo) 9 (21) 
1 271 

Having obtained $,(x,,x,) value in the coordinates of the controlled process, we can determine the value of 


$;(x,,X,) =1,535x, — 0,66 





Y; macrovariable in the coordinates of the system state and solve the third equation of the system (14). 
dx, db;(%,%) 
A agp ts Os Game) = 0. (22) 
After substituting the previously obtained 6,(x,,x,) values from (21) into (22) and using the system (13) at 


this step of decomposing the state space of the control system, we define the control value in the coordinates of the 
controlled process: 
PEE Th) GG ti): 2. 

= = 


TTT, © TTT, 
T,+T,+T, - 
25 066¢, ote) ggg BO 
DT, UD, 
The expression (23) determines the asymptotically stable control action on the control system described by 


U =1,5x, + 2x, +1,535x, +1,5x, 








(23) 





the equations (13). 
The simulation results of the system of differential equations (13) with the required value of the engine rotor 





angle: x;9=50 rad, and the values of 7, = 0.1,7, =0.2, 7, =0.3 time constants entered are shown in Fig. 4. 


Xi, rad 




















b) 





Fig. 4. Simulation results with control at 7, = 0.1, 7, =0.2 , 7, =0.3 : transition process on output coordinate (a); 


control change schedule (b) 


Fig. 4 shows that, in comparison with Fig. 3, the control structure is the same in both cases, that is, the point 
is that at the beginning, a positive control is applied, and then its sign changes to the opposite one. The settling time 
increases dramatically to the value close to 2.5 seconds, but with this, the maximum control reaches values greater 
than 5,000 V. This is unacceptable according to the maximum method. To limit the maximum control value, we 
introduce a nonlinear link of the “saturation” type with the same threshold as in the case of the maximum into the 
control loop. The results of the system simulation, with this threshold and with selected values of 7, 7,, T, 


constants, are presented in Fig. 5. 
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1 2 3 4 5 6 7 8 9 10 t. Sec 














Fig. 5. Simulation results with control at J, = 0.15, 7, =0.16,7, =0.17: 
output coordinate transition (a); control change schedule (b) 
Fig. 5 shows that the introduction of a nonlinear constraint on the control of the “saturation” type enables to obtain the 
desired result from the point of view of the control admissibility. In this case, the control structure determined by the 
maximum method is really observed, but the control is not time optimal. It should be noted here that such a 
modification of the ADAR method could also lead to the loss in system stability (Fig. 6). 


Xi, rad 



































b) 


Fig. 6. Simulation results with control at 7, = 0.015 , 7, = 0.016, 7, =0.017: 


output coordinate transition (a); control change schedule (b) 


However, despite the result presented in Fig. 6, the addition of amplitude constraint of the control signals to the 
control system obtained by the ADAR method essentially brings it closer to the solutions obtained through 
implementing the synthesis procedure by the maximum method. 

Discussion and Conclusions. Thus, for the basic model of the control system (2), the time optimal control was 
synthesized by the maximum method. For the same case, in the third part of the paper, with the same system parameters 
and boundary conditions, the control of the ADAR method was synthesized. Both of these controls look different. 
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In case of the maximum method, this is a software piecewise-constant control where the optimal time instant of 
the sign change of the control signal is important. In ADAR case, it is a continuous smooth functional control 
dependence on the coordinates of the system state. However, their implementation shows the structural similarity 
expressed in changing the control sign. After the ADAR method modification, the structural similarity has become even 
more visible. Note that the parametric similarity of the implemented controls can be achieved through selecting the time 
constants introduced by the ADAR method during the synthesis. This allows the authors to formulate the following 
hypothesis. 

- Selection of the values of 7,,7,,7, time constants introduced by the ADAR method in the synthesis 


procedure, in case of modification of the obtained control by the method proposed in the article, makes it possible to get 
time optimality of the synergetic control. The authors failed to obtain a general proof from this intuitive hypothesis. 
However, this task was not set within the framework of this paper. The numerical experiments with 7,,7,,7, variation 


helped us to achieve quasioptimality of the control obtained by the ADAR method. 
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